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Qh! Abstract 

We consider the CP violating asymmetries produced in the decay of heavy 

• i 

particles, studying the effects of heavy particle mixing for arbitrary mass splittings. 
■ A considerable enhancement of the asymmetries is achieved when the masses of the 

mixed states are comparable, and the enhancement is maximal for mass splittings 
of the order of the widths of the decaying particles. We apply the results to the 
particular case of heavy scalar neutrino decays relevant for leptogenesis scenarios. 



The violation of CP is one of the crucial ingredients identified by Sakharov [T|] as 
necessary conditions for a dynamical generation of the baryon asymmetry of the Uni- 
verse. Since the known CP violation in the standard model is probably too small to be 
helpful in generating the observed baryon asymmetry, the existence of new CP violating 
interactions in most of the extensions of the standard model is particularly welcome from 
this point of view. In the usual scenarios of baryon (or lepton) number generation by the 
out of equilibrium B (or L) violating decays of heavy particles, the CP violation arises 
in general from the interference of the tree level and one loop decay amplitudes, which 
allow the phases in the complex couplings involved to show up in the partial decay rate 
asymmetries. The one loop contribution which is usually taken into account is the so 
called 'vertex' one, in which two light particles produced in the decay of a heavy one 
exchange another heavy particle to produce the required final state (Fig. lb). However, 
many of the scenarios studied involve (and even require^) more than one 'flavour' of 
heavy particles, allowing then for a further possibility to produce CP violation. 

This new possibility arises from the so called 'wave-function' contribution |3|, |], |5], |5], 
0, in which a loop of light particles just mixes the initial state $ a with another different 
heavy state which then decays to the final state as shown in Fig. lc. This wave 
contribution turns out to be comparable to the vertex one when the heavy states have 
large mass splittings, and may be significantly enhanced for nearly degenerate states. 

The asymmetry in a global quantum number iV (for instance B or L) produced in 
the decay of a pair made of particle <3> a and its antiparticle $ a is given by 

/ 

where the quantity of interest to us is the partial rate asymmetry per decay into final 

1 The existence of just one heavy triplet scalar in the minimal SU(5) theory was actually considered a 
problem as regards the baryon number generation, since the CP violation appeared only at three loops 
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state / (of global charge Nf), given by 
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where BR($ a — > /) is just the branching ratio for the decay of particle $ a into the final 
state /. The wave function contribution to this quantity behaves, in the limit of large 
mass splittings, as ej a oc [M 2 — M 2 )' 1 , due to the propagator of the intermediate state 
$6, and hence is expected to be enhanced in the limit Mb — > M a . A general approach to 
study this quantity for arbitrary mass splittings has been considered in ref. ||. It is our 
purpose here to extend the formalism developed in ref. || and apply the results to the 
study of leptogenesis scenarios, for which a computation of the CP violation for large 
mass splittings was recently obtained JF]. 

To be specific, we will consider the case in which the heavy decaying particles are 
scalars, and ignore the vertex CP violation effects, which can be studied separately. The 
wave function mixing will have the effect of inducing an absorptive part in the heavy 
particle propagators, which will be responsible for the generation of the asymmetry. The 
effect of the one-loop self-energy diagrams in the propagators will be to modify the 
squared scalar mass matrix as follows 



where the renormalised mass matrix M 2 includes the dispersive part of the loops while the 
matrix T 2 arises from the absorptive part alone. The matrices M 2 and T 2 are hermitian, 
but H 2 is not. Hence, H 2 will be diagonalised in general by a non-unitary transformation 
matrix V [|J, i.e. 



This matrix V will then transform the initial 'flavour' states |$ a ) into the 'propagation' 
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eigenstates^ | i.e. 



& c ) = v- l \$ a ). 



(5) 



Similarly, for the antiparticle states |$ a ) ; one will have 



(6) 



These propagation eigenstates are the ones that will evolve as 
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Considering then the transition amplitude from the state |$ a ) to a final state |/), we 
have 



where Hi nt describes the interactions of $ a with the final state particles. From the 
superposition principle and using Eqs. (|3|,[§), one has 



The differential partial decay rate asymmetries arising from particle mixing will be pro- 
portional to the quantities 



It is interesting to notice that, in the limit of degenerate propagation eigenstates, i.e. 
uj c = u, these asymmetries vanish, as can be seen from eqs. (|9|,p!0|). 

To continue we will concentrate in the case of mixing between just two particles, for 

which the matrix V can be parameterised in terms of two complex mixing angles, 9 and 

2 The appearance of V~ x in eq. (Q) ensures that the kinetic term remains canonical, but the fact that 
V- 1 i^V^ implies that the propagation eigenstates are not orthonormal. 



T fa = (f\H int \$ a ), 




T fa (t)= n.cTfbV^Vcae 
f fa (t)= E^T/^y-V 



(9) 



A fa (t) = \T fa (t)f-\f fa (t)\ 2 - 



(10) 
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>, as follows 



rr_f cos9 —sin9e l ^\ 

V sinfle"^ cosfl ' ' ' ' 



Replacing this in Eq. (|) it is easy to obtain 

tt2 a tj2 tt2 

2i<t> n l2 . ( i . n .r,n\2 _ ^ n \2 n 2\ ( -, r>\ 

e =hi ■ {tg2B) -(hi-h:^' (12) 

where we recall that H 21 = M\ 2 — ^12- The eigenvalues of H 2 are then 

^1,2 = \ {#n + #22 ± ^Hli-HhY + ^Hli} ■ (13) 
After an explicit computation we then get 

A fa (t) = 2Re {T fl T* f2 [U la U* 2a - U a2 U* al }} + \T fb \ 2 {\U ba \ 2 - \U ab \ 2 } , (14) 

where we have defined 

U ab = V~ l W c V cbl (15) 

with 

W c ee e - ***. (16) 

In the case in which the initial state under consideration is an eigenstate of the (renor- 
malised) mass matrix, i.e. for M 2 b = M 2 5 ab , these expressions simplify considerably, and 
we have 

im(T n r; 2 r? 2 ) 

A fl (t) = 4 , 1 / /2 2|2 12J Re{K 2 - u 2 2 )(W; - W^XcosW! + sin 2 6W 2 )| , (17) 
— oj 2 I 

and a similar result holds for Aj 2 with the substitution W 2 «-> W\. 

We will then compute in detail the integrated rate asymmetry in this case. The 
branching ratios entering in eq. (0) are just 



//*oo 
dn a j Q dt\T fa (t)\' 



with dVL a the phase space element of particle $ a . We have then 

e% = N f tt a dtA fa (t), (19) 

J 

where fl a = M a /16n in our case of two body scalar decay. 
Integrating over time we find 

e w fa = 2N f n a lm {T fl T* f2 Y\ 2 } F a) (20) 

with 

{Ml -Ml) 



F 1 = Re{ul - ul) 

l^i ~ u 2 \ I 
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where we defined c<j a = m a — ?7 a /2. The result for F 2 is similar with the replacements 
7i ^ 72 and m\ <-> m 2 in the expression for 

In the case M| > Ml > |T 2 b |, one has 

2 

Fa ~> 7a (M|-M 1 2 )' (22) 
where in this limit 7 a ~ T 2 aa /M a becomes just the total width of particle $ a . Hence, the 
results obtained in the literature in the limit of large mass splittings can be recovered. 

For decreasing mass splittings, the function F x reaches a maximum, which for |T 2 2 | <C 
|r| 2 — r 2 x | takes place for 

Ml- Ml ~r 2 n + rl 2 . (23) 

The value of F 1 at the maximum is F t ~ Mx/^iY 2 ^ + Y 2 22 )}. 

On the other hand, for |r 2 x — Y 22 \ 3> |M 2 — M||, |r 2 2 |, i.e. in the limit of small mass 
splittings and small mixing, we get 

F ,, (Ml-M 2 1 ) 2M l \ AY 2 U Y 2 22 } 1 

1_ irfx-r 2 ,! 2 r 2 , \ (r 2 1 + r 2 2 ) 2 / M^ + ry [) 



This result coincides with the one in ref. || only in the limit <C T\ 2 (or T\ Y ^> r| 2 ) 
and if we neglect the last term, which although small is non-vanishing and survives in 
the limit Mf — > Mf (in which the propagation eigenstates are not degenerate due to 

^ r| 2 ). 

Another simple example is for the case Y\ x = r 22 , for which we have 



u\ - ujI = V(M| - Mf) 2 - 4|r? 2 | 2 (25) 

We see that here for |Mf — | = 2|r 2 2 | the two propagation eigenstates become degen- 
erate, and in fact for this mass splitting ej a = 0, but it will be non- vanishing otherwise. 
In particular, in the limitQ Mf -» Mf one has F x ~ -l/(2M 1 rf 1 ). 



Let us also emphasise that a crucial ingredient in all this computation is the proper 
specification of the initial state. The asymmetry of course depends on the starting ba- 
sis for $ a considered, and hence on the process which produces the initial state, so 
that ignoring the mixing at production would lead to incorrect results. For instance, 
in the case in which M 2 oc 1 (or more generally whenever the matrices M 2 and T 2 
commute), it is possible to change basis with a unitary transformation to make M 2 and 
T 2 both diagonal. Hence, the asymmetry computed in this new basis will vanish, since 
T 2 2 = now. However, these new states may not be the quantum states generated in 
the production process, and therefore the new basis may not be the appropriate one to 
compute the resulting asymmetry. 

Let us now consider the particular example of lepton number generation in the out 
of equilibrium decay of heavy scalar neutrinos, i.e. the supersymmetric version of the 



Fukugita and Yanagida scenario Ul0 |. In these type of models [IT, 0, [13], [13, II|, the 

3 The 'degenerate' situation rf x = Y\ 2 and M\ — M| actually is present in well known cases such 
as K°K° or B°B° mixing, where those constraints are imposed by CPT relations, and for which the 
integrated CP violating asymmetries are non- vanishing ||. 
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decay of the electroweak singlet (s)neutrinos, with masses M ^> TeV, produces a lepton 
asymmetry. This is then partially converted into a baryon asymmetry |16| by the effects 
of the anomalous B + L violation in the standard model [IT]], which is in equilibrium at 



temperatures larger than the electroweak phase transition one (~ 10 2 GeV). The study 
of the CP violation in these models, considering both the 'vertex' part as well as the 
'wave' contribution in the limit of large mass splittings, was carried out recently f7|. We 
now consider the effects of mixing for arbitrary masses using the formalism introduced 
above[]. 

The Lagrangian for the scalar neutrinos is, in a basis in which the mass matrix is 
diagonal, 

C = -\ ia e aP ^M a N* a LfH^ + (^P^fjV } + h.c. (26) 

where £f = (z/j ll) and H T = (H + H°) are the lepton and Higgs doublets (z = e, /i, r, 
and e a p = —ep a , with e 12 = +1). 

Since we are interested in the possible implications of small mass splittings, we will 
assume that the right handed neutrino masses consist of two almost degenerate states, 
with the third one being much heavier and hence effectively decoupled from the mixing 
mechanism. In this case, the effects of the third scalar neutrino can be included indepen- 
dently, and the mixing effects can be studied with the two flavour formalism discussed 
before. It is particularly interesting that scenarios with this type of spectrum have been 



widely considered in the literature [fL9|1 , and can naturally arise in 50(10) models. 



We will assume that sneutrinos are produced out of equilibrium by a certain unspec- 
ified mechanism (e.g. if sneutrinos are inflaton decay products [0 or the inflaton itself 

4 There has been a recent attempt to study the asymmetry in heavy Majorana neutrino decays in the 
limit of small mass splittings ]T^| , but those results are however at variance with ours. For instance, we 
find a dependence on Hyi through H12H21, and not through (H12 + H21) 2 as in ref. fL8| , although we 
expect similar results for neutrino and sneutrino decays in the supersymmetric model considered here. 
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[|I4}|), and for simplicity consider that the states produced initially correspond to one of 
the eigenstatesQ, say N\, of the mass matrix (so that M\ 2 = 0). Hence, the asymmetry 
will be given by eqs. fl2"0| , |2"l"| ) , where in this model a direct computation of the absorptive 



part of the sneutrino propagator leads to 

1 



r2 

ab 8tt 



^\) ba M a M b +(^\) ab s], (27) 



where the square of the four momentum will just be s = p 2 = M 2 by the on-shell 
condition. The first contribution to the r.h.s. of eq. ( p7| ) is given by the slepton and 
Higgs loop, while the second by the lepton and Higgsino one. 

As final states, we need to consider two possibilities, i.e. / = LfH 13 as well as 
/ = £"h cl3 . For the final state with sleptons, we have Tf a = —ie a ^\* ai M a l^fs, so that only 
the second term in the r.h.s. of eq. ( p7| ) contributes to the total asymmetry, and we have 

E i/I- {TnT}A} = ^Im {(AtA)L} = -^Im {(AU)^} , (28) 

where Lf = +1 is the lepton number of the final state. On the other hand, for the decay 
N a — > Ifh ^, one has Tf a = —ie a p\ ai) and only the first term in the r.h.s. of eq. ( pT|) 
contributes to the total asymmetry. Since now Lf = —1, we end up with the same 
contribution as the one coming from the slepton channel. (Notice that the asymmetry 
in a given final state channel results from the mixing generated by a loop involving the 
particles of the other final state channel.) 

So, summing the contributions from both final states we get 

^=-w Im { (A,A)ii K < 29 > 

where F\ is given in eq. (pl|). If we use that = (A^A)nMf /4tt, and the asymptotic 



5 If both Ni and N2 are simultaneously, but incoherently, produced, one needs just to add the asym- 
metries from both decays. 



expressions discussed previously for Mf ^> Mf 3> |r^|, one can see that in this limit 

£W = 1 M lM2 Im {( AtA )L} 

1 1-kMI-MI (AtA) n ' { } 

which coincides with the expression obtained in ref. J7|. 

In figure 2 we plot the total asymmetry for arbitrary mass splittings, normalised 
to the vertex contribution arising from the exchange of the second state N 2 J7| 

„ 1 /M|\ Im {( AtA )L} 

^ = -^ 9 {m!) (AtA) n ' (31) 

where g(x) = \fx ln[(l + 

Notice that e\ contains actually the same combination of Yukawas appearing in the 
wave function contribution ef (i.e. the factor ImKA^A)^}), so that the CP violating 
phase cancels in the ratio. 

In fig. 2 we adopted for definiteness T 2 22 = Mf/10, T 2 U = \T 2 12 \ = plotting 
the result as a function of x = M2/M1. In the limit of large mass splittings, the wave 
contribution approaches twice the value of the vertex one, as expected 0. For decreasing 
mass splittings (x — *■ 1), the enhancement in the wave contribution due to the mixing of 
the states is apparent, and reaches a maximum value ~ Mf/(r| 2 ln2) for Mf — M\ ~ Y\ 2 
(corresponding to x ~ 1.05 in this case) as discussed in eq. (|23|) . 

The dotted line corresponds to the asymptotic expression for the wave contribution 
in eq. (|3"0|), and gives a reasonable approximation to the result for Mf — Mf > 4X22 • For 
smaller values of Y\ % jM\ , the enhancement in the wave contribution is larger (and can be 
in principle of many orders of magnitude), and the asymptotic expression is valid down 
to smaller values of x. On the other hand, for Mf — M\ < r% 2 , £\ decreases significantly, 
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and for M2 — > Mi one has 

e i 1 11 



(32) 



which is tiny. The results are quite insensitive to the actual values of r 2 ^ and |r 2 2 |, as 
long as this last stays much smaller than Y\ 2 , so that the mixing angle 6 is small. If 
|rf 2 | ~ T\ 2 , the maximum enhancement is somewhat smaller but the general behaviour 
remains similar. 

It is important to keep in mind that the vertex and wave contribution arising from 
the exchange of the heavier third scalar neutrino N% may however be larger than the one 
coming from the exchange of the second one N2, even taking into account the possible 
enhancements for small mass splittings, due to the probably larger Yukawa couplings 
involved and the unknown size of the CP violating phases appearing in both channels 
(for three families, there are actually three independent CP violating phases entering in 



the asymmetry [20 ) . 



In conclusion, we have considered in detail the integrated CP violating asymmetries 
arising from heavy particle mixing, and studied the effects that appear when the mass 
splittings are of the order of the particle widths. The large enhancements which can be 
achieved can be helpful to explain the observed baryon asymmetry of the Universe, as 
we have exemplified with the study of a scenario for leptogenesis. 
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Figure captions: 



Figure 1: Diagrams which interfere to produce the CP violation in the heavy particle 
decay. Fig. lb gives the so called vertex contribution while Fig. lc gives the wave function 
one. 

Figure 2: Wave function contribution to the CP asymmetry, normalised to the vertex 
one, as a function of M 2 /M 1 , assuming M\ 2 = and taking T\ 2 = Mf/10, = |r^ 2 | = 

r| 2 /io. 
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